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Abstract 

We present a study of the elastic and plastic behavior of nanoporous gold in compression, focusing on molecular dynamics 
simulation and inspecting experimental data for verification. Both approaches agree on an anomalously high elastic 
compliance in the early stages of deformation, along with a quasi immediate onset of plastic yielding even at the 
smallest load. Already before the first loading, the material undergoes spontaneous plastic deformation under the 
action of the capillary forces, requiring no external load. Plastic deformation under compressive load is accompanied by 
dislocation storage and dislocation interaction, along with strong strain hardening. Dislocation-starvation scenarios are 
not supported by our results. The stiffness increases during deformation, but never approaches the prediction by the 
relevant Gibson-Ashby scaling law. Microstructural disorder affects the plastic deformation behavior and surface excess 
elasticity might modify elastic response, yet we relate the anomalous compliance and the immediate yield onset to an 
atomistic origin: the large surface-induced prestress induces elastic shear that brings some regions in the material close 
to the shear instability of the generalized stacking fault energy curve. These regions are elastically highly compliant and 
plastically weak. 
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1. Introduction 

Nanoporous gold (NPG) made by dealloying has at¬ 
tracted much attention due to its prospective applications 
in, e.g., actuation CHS], catalysis 0 m , sensing 0 mg, 
and microfluidics m- The material represents a net¬ 
work of short fibers or ‘ligaments’ with dimensions that 
can be controlled between 5 nm and several micrometers 
mu]. Significantly, dealloying conserves the crystal lat¬ 
tice of the parent alloy (in the order of 10 — 100 /im), so 
that ligaments are crystallographically coherent over large 
distances mm- Representing an ensemble of intercon¬ 
nected ligaments, NPG offers an implementation of the 
mechanical properties of nanowires, including specifically 
their high strength, into a macroscopic material. The ma¬ 
terial is thus a model system for study of the mechanical 
properties and the mechanisms of plastic deformation at 
the nanoscale. 

Experiments, and specifically the tests of macroscopic 
samples to large compressive strain of Refs. min hi], 
have established characteristic features of NPG deforma¬ 
tion which are now briefly summarized and compared to 
recent simulation studies using molecular dynamics (MD) 
use)] and finite element modeling (FEM) [20]. 
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NPG, while brittle in tension m m no ns, can be 
highly deformable in compression. Early experiments demon¬ 
strated this for focused-ion-beam (FIB) cut micropillars 
ESEIj, but improved synthesis protocols now allow de¬ 
formation to large compressive strain for mm-sized sam¬ 
ples mmm that are free of artifacts from FIB cutting. 
Gibson-Ashby scaling equations [25] |26] have been used 
for estimating the local yield stress of ligaments based on 
the effective materials behavior of NPG [23] [27] . However, 
the results of different measurements disagree. Micro and 
nanoscale tests suggest that the yield stress of ligaments 
increases with decreasing size, approaching the theoretical 
shear strength at sizes of a few nanometers [23j HU I27U30] . 
The results for mm-sized samples confirm the trend, yet 
they reveal systematically lesser strength HU [22]. This 
discrepancy partly reflects inconsistencies in converting 
nanoindentation hardness to yield stress [l6l 122] , More 
importantly, however, the compression experiments testify 
to an extended elastic-plastic transition, with early yield¬ 
ing and with a strong initial increase of the flow stress 
[2D]. The yield strength is therefore poorly defined. MD 
studies have so far not reproduced this finding, indicat¬ 
ing instead a pronounced linear elastic regime m and 
compressive yield strength in agreement with the scaling 
equations [19]. 

Most studies find NPG more compliant than predicted 
by the relevant Gibson-Ashby law [2HI EH E2]. There are, 
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however, two reports of unusually high stiffness at small 
ligament size |33j [34] . The recent tests on mm-sized sam¬ 
ples show an unusually high initial compliance and a sig¬ 
nificant stiffening during the early stages of plastic defor¬ 
mation nuns]- This evolution in elastic response cannot 
be explained by Gibson-Ashby type scaling of the effective 
stiffness with density [20] . While microstructural disorder 
in the network of ligaments was found to favor high com¬ 
pliance m , the mechanisms behind the unusually high 
initial compliance and the strong stiffening during early- 
stage deformation of NPG remain yet to be identified. In 
fact, large discrepancies between the reported experimen¬ 
tal stiffness data call for an independent verification. 

Experiments show practically no transverse strain dur¬ 
ing compressive plastic deformation under uniaxial load 
m- This observation is reproduced by MD simulations 
and FEM [19, 20., with FEM linking the observation to 
microstructural disorder [20]. As a consequence of the ab¬ 
sence of transverse strain, uniaxial compression densities 
the network. Consistent with this densification and with 
the scaling laws, compression brings an unusually strong 
strain hardening [T6] . It has been suggested that this im¬ 
plies a compression-tension asymmetry, with strain hard¬ 
ening in compression and strain softening in tension m. 
The tensile brittleness of NPG might thus reflect the plas¬ 
tic instability that ensues from strain softening. Indeed, 
constraining the volume change by infiltrating the pores 
of NPG with a polymer, while strengthening the material, 
suppresses the strain hardening and establishes ductility in 
tension El- Tensile deformation of NPG has been probed 
by atomistic computer simulation, with one study con¬ 
firming the tensile strain softening [18] while another finds 
strain hardening jT9]. The origin of this different behavior 
has not been addressed. 

Dislocations are accumulated during plastic deforma¬ 
tion of NPG and eventually even form cell walls. Electron 
backscatter diffraction detects the ensuing mosaic struc¬ 
ture at length scales larger than the ligament size [16] . 
while high-resolution transmission electron microscopy shows 
stacking faults and twins generated upon compression ffl. 
Strain hardening, besides partly reflecting densification, is 
also affected by interactions between the stored disloca¬ 
tions. This is revealed by the comparison of experiment to 
results from finite element modeling [20 and is also sup¬ 
ported by the emergence of strain-rate sensitivity as the 
defect density increases during deformation m ■ 

In this study, we present a study of NPG compres¬ 
sion deformation by MD simulations, addressing specifi¬ 
cally the immediate plasticity and subsequent strain hard¬ 
ening, as well as the unexpectedly high initial compliance 
and subsequent stiffening. For verification, experimental 
data are shown alongside the simulations. 


2. Methods 

2.1. Simulation 

A porous microstructure on a rigid crystal lattice was 
generated by spinodal decomposition of a binary mixture, 
using the Metropolis Monte Carlo algorithm with an Ising 
type Hamiltonian and periodic boundary conditions in all 
directions. A cubic box of 100 x 100 x 100 lattice spacings 
(408 A) and {100}-oriented edges was randomly filled with 
atom types A and B at 30 % atomic fraction A. Interac¬ 
tion energy for nearest neighbor atoms was zero for pairs of 
the same type, and 0.02 eV for pairs of different types. Af¬ 
ter the phase separation, all H-type atoms were removed. 
All isolated A-type atoms were also removed to assure a 
connected structure. As a result, the solid (volume- and 
mass-) fraction, <^, was 0.297. The same procedure was 
applied to H-type atoms to ensure open porosity. 

For the molecular dynamics simulations we used the 
open-source code LAMMPS [36] with an EAM potential 
for gold m- In most simulations the temperature was 
maintained at 300 K, while selected runs explored low tem¬ 
perature conditions (T = 0.01 K) to assess the conceivable 
role of thermally activated processes. 

In each run, the energy of initial structure was first 
minimized by the conjugate gradient method. To ensure 
a force-free structure after energy relaxation, the mini¬ 
mization was performed until the specified force tolerance 
is less than 10 -4 eV/A, while the energy criterion was 
set to zero. At the point where the minimization con¬ 
verged, the relative change in energy was less than 10 -12 . 
Then, the structure was thermally relaxed for 1 ns at zero 
load. Uni-axial compression was attained by scaling the 
whole simulation box in one dimension at each time step 
at the engineering strain rate of 10 8 s -1 , while keeping 
other dimensions at zero stress. Temperature and pres¬ 
sure were maintained by a Nose-Hoover thermostat and 
barostat [38 ? 39] with Martyna-Klein-Tuckerman modifi¬ 
cations @0741. The time step for all simulations was 2fs. 
All strains in this work are specified as compressive engi¬ 
neering strain, e. 

The evolution of Young’s modulus of NPG was ob¬ 
tained by carrying out another simulation, using simi¬ 
lar conditions but with intermediate unload/reload cycles. 
The effective Young’s modulus was determined as a tan¬ 
gent modulus m from the unload stress-strain segments. 

The open-source software Ovito m was used for vi¬ 
sualization. Dislocation defects in the material were de¬ 
tected and identified with the dislocation extraction algo¬ 
rithm (DXA) |44j 45 . This computational tool allows us 
to measure dislocation densities (broken down by Burgers 
vector) during the simulation. A surface reconstruction 
algorithm [46] was used to generate a geometric represen¬ 
tation of the evolving crystal surfaces from the atomic po¬ 
sitions and measure surface area and solid volume. Pla¬ 
nar defects such as stacking faults and twin boundaries 
were detected and classified by means of an atomic pat- 
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tern matching algorithm [46 that is based on the adaptive 
common neighbor analysis (a-CNA) method [47] . 

2.2. Experiment 

NPG specimens for the experimental study were syn¬ 
thesized by electrochemical dealloying as described in Ref. m 
The processing of the Au 2 sAg 75 alloy precursor for the 
dealloying included arc melting, homogenization at 900 °C 
for 120 h in vacuum, wire drawing, and cutting to form 
cylindrical samples of 1.19 — 1.21 mm in diameter and 
1.90 — 2.00 mm in length. These samples were then an¬ 
nealed in inert atmosphere (800 °C, 2 h) and used for dis¬ 
solution at 0.75 V in 1 M HCIO 4 electrolyte at room tem¬ 
perature. Following this, a polarization at 1.1 V until the 
current decayed to 10 pA and potential cycling at a scan 
rate of 5 mV/s between —0.5 V and 1.0 V for 20 cycles was 
performed on as-dealloyed NPG in fresh 1 M HCIO 4 solu¬ 
tion. This served to remove residual Ag as well as strongly 
bound oxygen from the surface. Finally, the NPG samples 
were immersed in deionized ultra-pure water for at least 
12 h and left for drying in Ar flow for more than two days. 
Scanning electron microscopy (SEM) yielded a ligament 
size of 40 ±5 nm in the as-prepared specimens. The initial 
solid fraction was estimated at ip = 0.298 db 0.01, based 
on the sample dimensions and mass. Similar to previous 
study Pi, the experimental sample is polycrystalline with 
no apparent crystallographic texture. With grain sizes in 
the order of 50 /im, the individual crystallites are three 
orders of magnitude larger than the ligaments. 

Compression testing used a mechanical testing rig (Zwick 
Z010 TN) at constant engineering strain rate 10 - 4 s -1 and 
under ambient conditions. A stress-strain curve was ac¬ 
companied by unload/reload cycles to obtain the effective 
stiffness. The minimal stress for each unloading step was 
1.75 MPa. A detailed description of the stiffness evalua¬ 
tion can be found in Ref. [SO]. 

The removal of oxygen species from the surface of NPG 
concurs with coarsening. As adsorbed oxygen has a strong 
influence on the mechanical behavior of NPG pang , we 
here compare mechanical response of the virtual sample 
with that of samples which, albeit having much bigger lig¬ 
ament size, exhibit clean surfaces. 

3. Results 

3.1. Morphology and initial relaxation 

Figure shows a snapshot of the virtual sample as 
generated by the Monte Carlo simulation (left), as well 
as its surface reconstruction (right). Similar to previous 
MD studies, where a phase-field spinodal decomposition 
approach was used for generating the boundaries of the 
ligaments mm, we find that our sample reproduces the 
interconnected nanoscale network structure of nanoporous 
gold well (see the SEM image in Fig. [1J). 

The specific surface area (area per volume of solid) of 
our virtual sample is a = 1.05/nm. In order to relate 
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Figure 1: (Color online) Nanoporous gold (NPG) structures used in 
this study, (a) As-prepared structure from spinodal decomposition 
(left) and its surface reconstruction (right). This sample has a rela¬ 
tive density of 0.297 and ligament size of 3.15 nm. Capillary forces 
induce plasticity, letting the sample volume decreased by 3.6 % after 
equilibration. Thus, lattice defects can be observed in the relaxed 
structure before deformation. Some examples of those defects are 
shown in (b). Throughout this work, if not otherwise stated, fee, 
twin boundary, stacking fault, and unclassified atoms are coded in 
blue, red, yellow, and green, respectively (see legend), (c) Scanning 
electron micrograph of the microstructure of a NPG sample. Rel¬ 
ative density is 0.298 =b 0.01, ligament size is 40 =b 5 nm. Note the 
similarity to the microstructure of the virtual sample. 


a and cp to a characteristic structural length-scale, L, we 
used the conversion rule: 

_ 1.63(1.25 - <p)(1.89 + <p(0.505 + ip)) 

L ~ a ’ (1) 

which yields L = 3.15 nm. This rule results from modeling 
NPG as a periodic diamond-like structure with cylindri¬ 
cal ligaments of diameter L, connected in spherical nodes 
which are just large enough to just cover the overlapping 
cylinder segments [20] . For the experiment, electrochem¬ 
ical capacity measurements [48] yield the volume-specific 
surface area of the tested sample, a = 0.0923 =b 0.002/nm. 
The conversion rule given in Eq. 0 thus puts the liga¬ 
ment diameter of the respective diamond-like structure at 
L = 35.8 zb 0.8 nm, which is well compatible with the value 
estimated from the SEM image. 

At ~ 40 nm edge length, our simulation box dimensions 
are one order of magnitude larger than the ligament size 
of ^ 3 nm. This suggests that the box size is sufficient 
to rule out finite-size effects in the MD simulation. For 
verifying that the simulation box contains a statistically 
representative volume of the virtual nanoporous structure 
generated by spinodal decomposition, we created several 
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virtual samples of same density and structure size by in¬ 
dependent runs of spinodal decomposition under identical 
conditions. All virtual samples were found to exhibit con¬ 
sistent mechanical behavior, supporting the adequate size 
of the box. 

Relaxing the initially rigid and defect-free porous crys¬ 
tal at T = 0.01 K simulation and with no external load led 
to 2.96 % decrease in the total volume, increasing the ef¬ 
fective solid fraction to ip = 0.306. Not a single dislocation 
could be found in the relaxed structure, but the relaxation 
generated a very small number of stacking faults and twin 
boundaries. The fault density was ^ 7 x 10 _5 nm _1 , corre¬ 
sponding to only ^120 atoms in faulted configurations, a 
negligible fraction of the ~ 1.2 million atoms in the entire 
structure. This suggests that the volume contraction after 
relaxation at low temperature is almost entirely elastic and 
can be understood as a natural consequence of the action 
of the surface stress. While the average surface-induced 
prestress can be calculated by using generalized capillarity 
equation (Section 4.1 below), the influence of surface stress 
is readily observable by monitoring the stress distribution. 
As an example, Fig. [2] shows a slice cutting through NPG 
sample before (left) and after (right) the initial energy 
minimization, with atoms color-coded according to the lo¬ 
cal von Mises stress (assuming equilibrium atom volume 
for all atoms). It is evident that although lattice defects 
have not shown up yet, the bulk atoms already experience 
changes in local stress, with accentuating changes near the 
surface or at the smaller segments of ligaments, implying 
that, before the onset of straining, the relaxed structure 
already suffers large stress (with non-uniform distribution) 
induced by the surface stress. We note that Fig. [2] agrees 
with previous MD study of NPG [19] where relaxation of 
initially unstrained NPG samples leads to changes in the 
corresponding histogram of local stresses. 
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Figure 2: (Color online) Surface stress influence on local stress 

distribution. This figure shows a slice cutting through NPG sample 
before (left) and after (right) the initial energy minimization, with 
atoms color-coded according to the local von Mises stress. Note color 
change in the bulk regions upon energy relaxation, indicating local 
shear stress balancing the surface stress. 

As compared to the relaxation at T = 0.01 K, equili¬ 
brating the starting structure at 300 K brings more pro¬ 


nounced structural changes. Figure ID> illustrates that 
many lattice defects are generated; the densities of dis¬ 
locations, stacking faults and twins are 4.8 x 10 -4 nm -2 , 
3.2 x 10 -3 /nm -1 and 1.4 x 10 -3 nm -1 , respectively. The 
creation of lattice defects even in the absence of an ex¬ 
ternal load is also observed in experiments. Transmission 
electron microscopy (TEM) reveals this process for NPG 
when ligaments yield to the surface-induced stress during 
fast dealloying m- 

The 300 K equilibration increases the solid fraction to 
ip = 0.308. The corresponding specific surface area and 
ligament diameter are a = 1.04/nm and L = 3.16 nm, 
respectively. The equilibration creates a stable structure, 
and in particular - apart form the generation of lattice 
defects - the initial long-range coherent crystal lattice is 
maintained. Bearing in mind the periodic boundary con¬ 
ditions, one can consider the virtual NPG structure as an 
infinitely extended single crystal with pores and disloca¬ 
tions. 

3.2. Stress-strain behavior under compression 

Next we studied stress-strain behavior in compression. 
As a background for inspecting the simulation results, Fig. 
shows the experimental stress-strain curve for a mm-sized 
NPG sample with L = 40±5 nm described above. The sim¬ 
ulated stress-strain curves of the virtual sample at 300 K 
are shown in Fig. Hf, with the continuous loading in red 
and the load/unload sequence in blue. For the purpose of 
clarity, data at strains <0.1 are magnified in the insets. 

The experimental data in Fig. [3^i are well compati¬ 
ble with data for a sample of similar L reported earlier 
in Ref. [16] (cf. Fig. Id there). Remarkably, despite the 
deviation in ligament size between experiment and sim¬ 
ulation, the respective stress-strain curves show excellent 
agreement. The small-strain deformation is nonlinear with 
a quasi immediate onset of plastic deformation already at 
the smallest load. This immediate yielding is suggested 
by the residual strain which is found upon unloading even 
after small deformation. 

Experiment and simulation also agree in relation to the 
unusually pronounced strain hardening at larger strain. 
For example of the virtual sample, the flow stress under 
continuous loading is as low as 18.2 MPa at strain 6 = 0.1, 
rising to 93.1 MPa at strain £ = 0.4 and finally attaining 
591 MPa at the end of compression. 

The stress-strain curves of the virtual sample at T = 
0.01 K are shown in Fig. [4] with graphs color-coded as in 
Fig- Ha- The stress-strain curves of virtual NPG are found 
for all practical purposes independent of the temperature. 

We emphasize the deformability of NPG, real and vir¬ 
tual, in compression, along with low yield strength and 
high strain hardening. To elucidate the underlying mecha¬ 
nisms, we now inspect the deformation-induced structural 
changes, concentrating on the deformation in continuous 
compression. 
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Figure 3: (Color online) Deformation curve of nanoporous gold, (a) 
Experimental compressive stress-strain curve of a millimeter-sized 
sample subjected to a sequence of loading-unloading, (b) Simulated 
stress-strain curves of virtual sample at 300 K. While the red curve 
depicts behavior of this sample in continuous compression to the axial 
engineering strain of 80 %, the blue one presents its corresponding 
behavior in a sequence of load/unload. The inset magnifies the load- 
unload curve at strains less than 10%. Throughout this work, cr and 
£ denote compressive stress and axial engineering strain, respectively. 


3.3. Transverse strain and surface area 

Figure [5^i shows the evolution of transverse strain dur¬ 
ing compression in our uniaxial load simulation. It is seen 
that the lateral dimensions remain essentially invariant for 
e up to 0.5, which is in agreement with experimental data 
measured during plastic flow [HI ITS], Consequently, the 
change in axial dimension translates directly into volume 
reduction. In other words, compression entails densifica- 
tion. By the end of straining (at 5 = 0.8) the solid fraction 
of the virtual sample attains the value of ip = 0.918. 

One of the signatures of microstructural evolution is 
the variation of the specific surface area. The graph of 
a - scaled to the initial value, <ao, of the structure before 
straining - is shown in Fig. [5]3 (entire evolution of recon¬ 
structed surface during deformation is shown in Movie-1 
in Supplementary data). The data set reveals a continuous 
loss of surface area during compression. 
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Figure 4: (Color online) Deformation curve of nanoporous gold. 
These simulated stress-strain curves were recorded at T — 0.01 K. 
As in Fig. |3j), the red and the blue curves represent the behavior 
under continuous loading and load/unload sequence, respectively. 
While the surrounding temperature differs significantly (0.01 K vs. 
300 K), similar behavior is obtained, as is clearly seen here and in 
Fig- [3 )p 


3.4. Evolution of effective Young’s modulus 

The load/unload segments of the experimental and sim¬ 
ulated test yield data for the effective Young’s modulus, 
Y eff , which we computed as an unloading tangent mod¬ 
ulus [42]. We verified that secant moduli, which average 
over the unload/reload process, provide consistent results. 
As the virtual samples density during compression, we 
found it of interest to plot Y eff versus the solid fraction. 
Figure [f^i shows Y eff and ip in log-linear plots, revealing 
that the virtual sample starts out exceedingly compliant 
at both temperatures of the MD simulations. The stiffness 
then increases substantially as the compression proceeds. 
For example, at room temperature, while the first value of 
Y eff at 5% strain is 542 MPa, the last value of Y eff at 57% 
strain is 9.03 GPa, nearly 20-fold stiffer. 

The experimental effective stiffness is shown along with 
the simulation data for comparison (see Fig.[6j circles). In 
view of the unusual compliance found by MD we empha¬ 
size the remarkable precision with which the simulation 
matches the experimental data, particularly at low strain 
where the compliance is highest. 

The relation between effective modulus and relative 
density of macroscopic metal foams is the subject of one 
of the Gibson-Ashby scaling laws [25]: 

yeff = ybulk^2 ^ ( 2 ) 

where Y bulk denotes Young’s modulus of the bulk metal. 
We emphasize that this equation describes the variation in 
stiffness with the initial density of the foam structure. It 
is well established that the stiffness variation with density 
during compression follows quite different trends, see for 
example Ref. [49]. Nonetheless, it is of interest to compare 
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Figure 5: (Color online) Structural change of NPG in compres¬ 
sion. (a) Change of transverse strain with respect to axial strain 
of nanoporous gold in compression. Lateral expansion is negligible 
for strain up to 50 %, as illustrated in the inset, (b) Change of sur¬ 
face area per solid volume, a, of nanoporous gold in compression. 
A continuous loss of surface area is observed, as shown in Fig. [5j). 
All values are scaled to the initial value of the undeformed structure, 
cko- Entire evolution of reconstructed surface during deformation is 
shown in Movie-1 in the Supplementary data. 

our data to Eq. © as a reference for the behavior of porous 
metals. As a basis we calculated the Young’s modulus of 
nontextured massive polycrystalline Au for the interatomic 
potential of our simulation at OK via Kroner’s formulation 
[50] . The result, Y bulk = 78 GPa, agrees well with the 
experimental value of massive polycrystalline Au at room 
temperature, Y bulk = 80 GPa [51 . Therefore, we take 
ybuik _ yg GPa and Y bulk = 80 GPa as the references for 
Y bulk at T = 0.01 K and room temperature, respectively. 

As an illustration for the first value of Y eff at 300 K, 
we substitute the corresponding solid fraction, = 0.321, 
into Eq. ([2|. The result is 8.2 GPa for Y eff , fifteen times 
higher than the value indicated by our simulation at that 
density. Consistently, at the first experimental value of 
the density, p = 0.299, the scaling law overestimates the 
experimental stiffness, for which we find 347MPa, by more 
than the factor twenty. 

Figure shows the evolution of Y eff normalized to 


Figure 6: (Color online) Change of effective Young’s modulus, Y eff , 
of NPG during compression. While Fig.[6ji shows the absolute values 
of Y eff versus relative density, cp, in a log-linear plot, scaled values of 
Y eff (to the Young’s modulus Y bulk of massive polycrystalline gold) 
are shown versus (p in a linear-linear plot (Fig. [6J)). Gibson-Ashby 
scaling law is also superimposed in Fig. [6]). The agreement between 
simulation and experiment is excellent. As is clearly seen in this plot, 
the evolution of Y eff is characterized by two features: exceptionally 
high initial compliance and ensuing stiffening under compression. At 
any time during compression, NPG is much more compliant than the 
prediction of Gibson-Ashby law. 

ybulk ^ i nc i uc q n g the graph of Eq. [ 2 ] superimposed to the 
simulation results. The figure illustrates that the effective 
stiffness of NPG is consistently very much lower than that 
of the scaling law. It is not unusual for metal foams to 
be more compliant than the prediction by the scaling law, 
see for instance Ref. [52], which reports a deviation by the 
factor three, and references therein. Yet, the magnitude 
of the deviation, which is apparent in our data (Fig. ©>) 
is extraordinarily high. 

In relation to the trends in Fig. If we also find it re¬ 
markable that the stiffness values of our real and virtual 
samples appear to extrapolate to zero at a finite solid frac¬ 
tion, slightly less than 0.3. This trend is clearly distin¬ 
guished from the power-law scaling embodied in Gibson- 
Ashby type laws, and it illustrates the extraordinary com¬ 
pliance of NPG. 
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3.5. Evolution of defect density 

Figure shows the evolution of dislocation density 
during compression. As is illustrated in the inset, Shock- 
ley partial dislocations - generated during relaxation as 
stated earlier in this work - are present already before the 
onset of loading. Compression then brings a continuous ac¬ 
cumulation of dislocations. As the deformation proceeds 
the fraction of full dislocations increases, yet the Shockley 
partials remain predominant throughout at all strains. 
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Figure 7: (Color online) Evolution of fault density, (a) Accumula¬ 
tion of dislocation density during compression at 300 K. Most con¬ 
tribution is from Shockley partials. (b) As a result of dislocation 
activity, twin boundaries and stacking faults populate during com¬ 
pression. The insets magnify parts of curves at strain less than 10 %. 
As is clearly seen from these insets, dislocations are already present 
before the onset of straining, and get their activity immediately upon 
loading. 


The largest dislocation density, at the end of the de¬ 
formation, is ^ 4 x 10 17 m -2 or ^ 0.4 nm“ 2 . Taking the 
characteristic volume of ligaments in the order of L 3 , we 
infer few (in the order of 4) dislocations per ligament. This 
number is well compatible with the snapshot of defect con¬ 
figuration on a cross-section through a virtual sample at 
55% engineering strain, Fig. [8j 

The inset in Fig. [7^ shows that the dislocation density 
increases immediately as the deformation sets in, even in 
the initial stages of compression. This observation rep¬ 


resents another manifestation of the immediate yielding, 
confirming the conjectures obtained from the absence of 
linear elastic stress-strain relation and from the lasting 
plastic deformation after unloading from few percent net 
elastic-plastic deformation. 

Figure [7 Jd is analogous to Fig. [7^i, except that the den¬ 
sities of twin boundaries and stacking faults are plotted. 
The analogous observations apply, with both densities in¬ 
creasing right from the start of loading and throughout 
the deformation. 

A typical defect structure under compression is shown 
in Fig. [8] (entire evolution of defect structures during de¬ 
formation is illustrated in Movie-2 and Movie-3 in Supple¬ 
mentary data). In line with previous MD studies [T8l [19] . 
many twin boundaries (red) and stacking faults (yellow) 
are observed. Lomer-Cottrell locks are also found, as il¬ 
lustrated in the lower magnification in Fig. [8] 
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Figure 8: (Color online) Typical defect structures found in 
nanoporous gold under compression. Left plot shows a slice cut¬ 
ting through sample at 55 % strain. Loading direction is perpendic¬ 
ular to plot plane. Magnifications of circled regions are shown on 
the right, depicting examples of grain boundaries (top) and Lomer- 
Cottrell locks (bottom) found in the deformed structure. Plasticity 
happens not only at the ligament junctions, but also in the middle 
of ligaments. Color coding is shown directly within the figure. En¬ 
tire evolution of defect structures during deformation is illustrated 
in Movie-2 and Movie-3 in the Supplementary data. 


Figure [5] also reveals that, in contrast to some experi¬ 
mental observations dans!, planar defects were found not 
only in the “nodes” which connect the ligaments, but also 
in the middle of the ligaments. 

At larger strains, grain boundaries form during densifi- 
cation, as exemplified by the upper magnification in Fig. [5] 
These boundaries result from cold-welding of neighboring 
ligaments. While this process contributes to densitying 
the structure, it is remarkable that the pore network, and 
with it the nanoporous structure, remains well-defined up 
to very large compressive strain, in spite of the largely 
increased solid fraction. 
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4. Discussion 

Let us start out the discussion of the stress-strain be¬ 
havior of NPG by a comparison to the behavior of macro¬ 
scopic metal foams. Experiment, for instance Ref.s m 
E 2 HM], on these materials suggests stress-strain curves in 
compression quite similar to NPG. Specifically, foams with 
higher initial density show the same type of progressive 
hardening with ongoing deformation as NPG. However, 
closer inspection reveals significant differences. Firstly, 
conventional metal foams deform elastically at small stress, 
contrary to NPG. Secondly, load-unload tests on conven¬ 
tional foams show the effective stiffness to remain roughly 
constant during compression or to even slightly decrease in 
its initial stages [49|. This contrasts with the pronounced 
stiffening of NPG upon compression that is observed in our 
study. Thirdly, compression tests of conventional metal 
foams of low initial density, comparable to our samples, 
exhibit a plateau in flow stress while the deformation is 
carried by localized crush bands [26l|52l [55]. NPG lacks 
the plateau, and in fact an experimental investigation of 
its microstructure evolution during compression testifies 
to uniform deformation as opposed to crush bands m- 

Our discussion of the distinctive behavior of NPG em¬ 
phasizes the role of its nanoscale microstructure. As one 
related aspect, it is well established that the strength of 
nanoscale objects typically increases with decreasing size 
f56ff58] and high strength comes along with low strain 
hardening [59 . As a network of nanoscale ligaments, NPG 
is expected to exhibit a similar trend in response to ex¬ 
ternal load. Yet, in line with previous experimental study 
m, our results show that both real and virtual NPG sam¬ 
ples are highly deformable, with ill-defined yield strength, 
extremely high compliance and pronounced strain harden¬ 
ing. Here, we address the most relevant issues. 

4-1. Surface-induced stress and initial dislocation content 

The microstructure of our virtual sample was generated 
by simulating spinodal decomposition on a rigid crystal 
lattice. By construction, the initial structure was there¬ 
fore free of strain and of lattice defects. Prior to the MD 
studies of deformation behavior, the virtual sample was 
subjected to a relaxation free of load. The ensuing den- 
sification by plastic deformation, accompanied by the cre¬ 
ation of various lattice defects, may be understood as the 
consequence of surface stress. The surface of a solid ex¬ 
erts forces in the underlying bulk (see below) that need to 
be compensated by bulk stress. For elongated structures, 
such as the ligaments of NPG, the surface-induced bulk 
stress has a shear component [60 that may prompt dislo¬ 
cation nucleation and spontaneous plasticity. Experiments 
on NPG detect the associated volume shrinkage and the 
generation of dislocations m, and atomistic simulation 
has previously confirmed the mechanism [BT|. 

While Fig. [2] documents the formation of these surface- 
induced stresses during the energy relaxation of the ini¬ 
tially unstrained NPG sample, the magnitude of the surface- 


induced stress, S, in the bulk of the nanoporous material 
is governed by the generalized capillary equation for solids 

m, 

SdF + /scL4 = 0, (3) 

where V, A , and s denote volume, surface area and surface 
stress tensor, respectively. This equation embodies a trend 
for even isotropic surface stress to create anisotropic stress 
in bulk, whenever the surface orientation distribution func¬ 
tion is anisotropic. Approximating the surface stress as 
isotropic with magnitude /, so that s = / P with P the 
surface projection tensor [62] . and considering a cylindri¬ 
cal nanowire as a cartoon of the ligament geometry, one 
obtains the principal values of surface-induced bulk stress 
in axial and radial directions, respectively, as Sa = —4 f /L 
and Sr = 2 f/L |60]. The maximum projected shear stress, 
of magnitude r = 2 f/L, is found on planes inclined by 
45° to the wire axis. For the EAM potential used in this 
study / ~ 1.1 J/m 2 [63] . With a ligament diameter of 
L = 3.15 nm, the surface-induced shear stress is thus esti¬ 
mated at r ~ 0.70 GPa. While this value is well below the 
theoretical shear strength of our potential, the disordered 
nature of the nanoporous network entails variations in lig¬ 
ament diameter and a fraction of local cross-sections that 
are well below 3.15 nm in diameter and that, therefore, 
may experience local surface-induced stresses of sufficient 
magnitude to induce spontaneous plastic shear. 

4.2. Early stages of plasticity under load 

The surface-induced plastic deformation during the equi¬ 
libration has important ramifications for the yield strength 
of the material. Most obviously, at room temperature, the 
virtual samples start out with a significant initial dislo¬ 
cation density. It has been well documented that the me¬ 
chanical behavior of nanoscale objects depends crucially on 
this parameter. While practically ideal strength may be 
achieved in dislocation-free whiskers, irrespective of their 
size, nanoscale pillars with an initial dislocation content 
are found considerably weaker [ 64H66] . The immediate 
yield onset of NPG may thus be promoted by its high ini¬ 
tial dislocation density. However, the very same behavior 
observed at nearly zero temperature where no dislocation 
was stored in the virtual samples indicates that the pre¬ 
existing dislocations are not the only reason for the imme¬ 
diate yielding. 

Equally relevant is the observation that not all liga¬ 
ments yield under the action of the surface stress dur¬ 
ing the initial equilibration. This is naturally understood 
as the consequence of a distribution of resolved surface- 
induced shear stresses. Only some ligaments will yield, 
thereby thicken and so reduce their local surface-induced 
prestress. Other ligaments may retain their prestress. It is 
then natural to assume that some ligaments are prestressed 
close to yielding, with only small external load required to 
locally reach the theoretical shear stress. In other words, 
the heterogeneous nature of the surface-induced prestress 
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contributes to early yielding and, when considering the 
entire porous network, to an extended elastic-plastic tran¬ 
sition as the external load increases. 

The above conjecture is perfectly born out by the ob¬ 
servations on the mechanical behavior of our virtual sam¬ 
ple and by the experimental data in Fig. [3] The last¬ 
ing plastic deformation after unload even from strain as 
small as ^ 1% testifies to the quasi instantaneous onset of 
yielding as load is applied, in agreement with the presence 
of very large surface-induced prestress in some ligaments. 
The same picture may also explain the initial increase in 
flow stress, which may reflect the exhaustion of local con¬ 
figurations with a low yield stress. In other words, as the 
deformation proceeds, the external load needs to increase 
in order to overcome those configurations in which the 
prestress is less active in assisting the deformation. Simi¬ 
lar behavior has been discussed in massive nanocrystalline 
metal when deformation events are confined to individual 
grains with a distribution of crystal orientation and hence 
a distribution of local Schmid factor m- Stress-strain 
curves of such nanocrystalline metal samples exhibit an ex¬ 
tended elastic plastic transition, yet with a preceding, well 
expressed elastic regime. The complete lack of a regime of 
purely elastic deformation is characteristic of nanoporous 
gold. 

The disordered nature of the porous network is another 
reason for stress heterogeneity in nanoporous gold. De¬ 
pending on its orientation relative to the macroscopic load, 
and depending on the load-bearing capacity of its neigh¬ 
boring ligaments, each ligament will suffer compressive 
load or bending moments which are related to the macro¬ 
scopic load by a local stress concentration factor. The wide 
distribution of stress concentration factors contributes to 
the gradual yield onset. This distribution should be well 
approximated by the FEM simulation of NPG in Ref. m- 
Yet, the simulated stress strain curves of that study do not 
reproduce the immediate yielding of the present work. The 
existence of an extended regime of purely elastic deforma¬ 
tion in the FEM simulation thus argues against distributed 
stress concentration factors as the origin of early yielding. 
These considerations again emphasize the importance of 
the surface-induced prestress that was discussed above. 

As another point of support for the above argument we 
advertise its implication of a compression-tension asymme¬ 
try. The surface-induced bulk stress would enhance com¬ 
pressive load but counteract tensile load. The strength 
should thus be enhanced in tension. This was pointed out 
by Farkas et al. Pi, who indeed confirmed the asymmetry 
in their MD simulation. In this relation we point out that 
a well-expressed linear elastic regime is seen in MD tension 
tests of nanoporous gold, but when compression is simu¬ 
lated with the same potential the elastic regime is largely 
diminished m- This is again compatible with surface- 
induced prestress as the origin of yielding in compression 
but not in tension. 

Surface roughness is another factor that might con¬ 
tribute to the deformability of NPG. For nanowires and 


nanopillars, MD simulations suggest that step edges and 
corners result in stress concentration that promotes dis¬ 
location nucleation and thus reduces the yield strength 
[68HZQ]. Recently, the local strain concentration at edges 
of nanostructures could be experimentally confirmed by 
imaging with aberration-corrected TEM m • In view of 
the impact of roughness for yielding it is significant that 
Fig- [D 3 reveals an abundance of step edges and corners in 
NPG; indeed at least part of these defects are geometri¬ 
cally necessary in order to achieve the nanoscale surface 
curvature of the ligament network. 

4-3. Strain hardening 

The comparison between the evolution of flow stress 
with strain in experiment and simulation (Fig. [3| reveals 
a remarkable agreement between both data sets. This sug¬ 
gests that the simulation catches the essentials of the strain 
hardening behavior of nanoporous gold. 

The strain hardening of NPG may partly be under¬ 
stood as a manifestation of densification under compres¬ 
sion. As pointed out in Section 3.3, the lack of transverse 
plastic strain implies that compressive axial strain entails 
a numerically equal volume compression and, hence, an en¬ 
hanced solid volume fraction. In view of the Gibson-Ashby 
scaling equation for the variation of flow stress with initial 
density [25], 

<7 flow oc <y/ 2 , (4) 

it is natural to also expect a qualitatively similar behavior 
when the material is densified during compression. In a 
constitutive sense, the increase in flow stress with ongoing 
strain represents an (effective, macroscopic) strain hard¬ 
ening. Of course, the link to the scaling law Eq. 0 can 
only be qualitative in view of the dissimilarities in struc¬ 
tures emanating from changing density during synthesis or 
during compression, respectively. In this qualitative sense, 
the expectation is supported by the observation that sup¬ 
pressing the volume change by infiltrating the pore space of 
NPG with a polymer will indeed eliminate the strain hard¬ 
ening El- Yet, the finite element simulation in Ref. f20j 
demonstrates the need to account additionally for work 
hardening in the constitutive equation of each nanoscale 
ligament, on top of the mere densification effect. 

Experiment points towards dislocation storage and dis¬ 
location interaction as the origin of the local work harden¬ 
ing within each ligament. Electron backscatter diffraction 
data of the microstructure evolution during deformation of 
nanoporous gold testify to the gradual evolution of mosaic 
spread at a scale larger than the ligament size [16] . This 
requires dislocation storage, and the very large increase 
in dislocation density in our simulation clearly supports 
this conjecture. What is more, the experimental data re¬ 
veal a gradual increase of the strain rate sensitivity as 
the deformation proceeds, an observation that is consis¬ 
tent with an increasing role of dislocation interaction m- 
Our simulation supports that picture as well. The final 
dislocation density of several dislocations per ligament is 
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extremely high, suggesting drastic Taylor hardening. The 
Lomer-Cottrell locks that are apparent in our deformed 
structure (as exemplified in Fig. [8| provide a direct evi¬ 
dence of dislocation interaction. We note that the contin¬ 
uous accumulation of lattice defects during compression 
and the impact of dislocation interaction implies that the 
dislocation-starvation scenario, which is often considered 
in models of small-scale plasticity mm, does not apply 
for NPG. 

The strain rate of MD simulation is much higher com¬ 
pared to that of our experiment. When there are com¬ 
peting deformation mechanisms with different strain-rate 
sensitivity, as in nanocrystalline metals mm, the dom¬ 
inant deformation mechanism depends on the strain rate. 
Specifically, classical MD simulations may overestimate 
the contributions of dislocation slip m • In that respect, 
we emphasize that our simulation results as well as pre¬ 
vious studies mmm consistently highlight dislocation 
slip as the only governing process in deformation of NPG. 
We therefore expect this mechanism to dominate indepen¬ 
dently of the strain rate, supporting a discussion of simu¬ 
lation and experiment within the same conceptual frame¬ 
work. 

4-4- Initial compliance and subsequent stiffening 

As is presented in Fig. [6j the evolution of Y eff is charac¬ 
terized by two features: exceptionally high initial compli¬ 
ance and a subsequent stiffening under compression. The 
excellent consistency between simulation and experiment 
indicates that the unusual behavior is indeed intrinsic for 
NPG. 

As with the early yielding, the elastic softening must 
partly be related to the microstructural heterogeneity. Due 
to the topological disorder, different regions will exhibit 
different elastic response, and soft regions might contribute 
to high effective compliance. FEM modeling [20] con¬ 
firms the strong correlation between structural disorder 
and compliance in NPG. Yet, the FEM model fails to re¬ 
produce the strong stiffening during the initial compres¬ 
sion. This points towards the role of atomic-scale pro¬ 
cesses which are not reflected in the constitutive laws of 
the FEM. 

A conceivable atomistic origin of the anomalous com¬ 
pliance is reversible microplasticity based on the bow-out 
of pre-existing dislocation segments between their anchor 
points. For an example of this effect, here in cold-worked 
iron, see m- Yet, the equally high compliance of the 
essentially defect-free virtual NPG of our low-temperature 
runs immediately rules out a decisive impact of preexisting 
dislocations on the compliant response of NPG. 

Similarly to the elasticity of nanowires, a more con¬ 
sistent explanation of the anomalous compliance consid¬ 
ers surface-related phenomena. Yet, the origin of size- 
dependent effective stiffness values of nanowires is contro¬ 
versial, see e.g. Ref. [78, and references therein. Surface 
excess elasticity is an obvious explanation US ES, and 
is supported by the observation of large changes in the 


stiffness when the surface state is reversibly varied [35] . 
yet an agreement on the sign of the surface excess elas¬ 
tic constants, stiffer or more compliant than bulk, has 
yet to emerge. The role of elastic nonlinearity of the 
core of wires prestressed by the capillary forces has also 
been advertised (e.g. see [8lH83] b the issue is challeng¬ 
ing since the nonlinear elastic response of cubic crystals 
is complex, see e.g. [S3]. Irrespective of their details, the 
above considerations require the modified stiffness to de¬ 
pend on the specific surface area. It is therefore remark¬ 
able that our findings show substantial stiffening of NPG 
after small plastic strain, while the net surface area is es¬ 
sentially invariant. This points towards processes other 
than surface excess elasticity or modified elastic response 
of the bulk. One issue stands out in the present geom¬ 
etry: as discussed in Section 4.1, the surface stress in¬ 
duces shear stresses which exceed the theoretical shear 
strength in some ligaments and closely approach it in oth¬ 
ers. Since the theoretical shear strength is determined by 
a point of inflection in the generalized stacking fault en¬ 
ergy curve, it entails a vanishing shear modulus. This is 
illustrated schematically in Fig. [9] In other words, the 
surface-induced stress and the concomitant shear will sys¬ 
tematically reduce the shear stiffness of the ligaments in 
NPG. This is well consistent with the anomalous com¬ 
pliance of the material in its initial state. Furthermore 
one expects that, as the material is plastically compressed, 
more and more of these near-unstable configurations will 
be pushed over to stable faulted or twinned configurations, 
relaxing their stress and reverting back to a configuration 
of conventional, high shear stiffness. This agrees with the 
observed stiffening during the early stages of plastic defor¬ 
mation. Even though the shear instability in this picture 
seems unusual, it is well supported by all observations of 
our work. 




Figure 9: (Color online) Schematic illustration of shear instability. 
Graph at left shows generalized stacking fault energy function, 7 , ver¬ 
sus shear displacement, u, between adjacent atomic planes. Graph at 
right shows the associated variation of the shear modulus, G, which 
scales with the second derivative of 7 (u). The negative-valued second 
derivative between the points of inflection of 7 (u) implies instability 
to shear. 


4.5. Size-dependent flow stress 

The discussion so far advertises a truly remarkable 
agreement between simulation and experiment, including 


10 




the stiffness and its evolution during compression, the ini¬ 
tial yield behavior and the strain hardening. This agree¬ 
ment supports the notion that the model structure is com¬ 
patible with experiment in respect to its topological fea¬ 
tures and likely also in respect to the state of the sur¬ 
face. Yet, the agreement breaks down when it comes to 
the structure size, where the ligaments of the experimental 
sample are around tenfold larger than those of the simu¬ 
lation. As a consequence, and in view of the well-known 
trend for strength enhancement at reduced size, the nu¬ 
merical agreement between the flow stress values in simu¬ 
lation and experiment is surprising. As the agreement be¬ 
tween the stiffness values suggests that the topology of the 
simulation is realistic, we suspect that the unexpectedly 
low flow stress of the simulation is linked to the surface 
geometry. Specifically, Fig. [TJd illustrates that the small 
structure size in the simulation entails extremely high local 
surface curvatures and, in the consequence, extreme den¬ 
sities of geometrically necessary step edges and kinks. We 
have already mentioned in Section 4.2 above that atom¬ 
istic simulation studies EHHZO] point towards a strong in¬ 
fluence of surface defects on the yield behavior. The low 
flow stress in the simulation supports these remarks, sug¬ 
gesting that surface roughness has an important impact 
on local stress concentrations and, hence, on the macro¬ 
scopic yield behavior of the extremely small structures of 
our simulation. 

5. Summary 

We have presented a study of the mechanical behav¬ 
ior of nanoporous gold (NPG) in compression by means of 
molecular dynamics simulations. A virtual sample pre¬ 
pared by spinodal decomposition with ligament size of 
3.15 nm and solid fraction ip = 0.297 was uni-axially loaded 
and thus strained to very large deformation. Experimental 
data are shown for comparison. 

Our study finds close agreement between experiment 
and simulation, advertising the following common phe¬ 
nomenology: 

• NPG can be deformed to large strain in compression, 
with very high strain-hardening coefficient. 

• NPG yields immediately even at the smallest applied 
load. 

• The plastic Poisson ratio remains essentially zero 
throughout a large range of strain. 

• As the most remarkable observation, NPG is anoma¬ 
lously compliant, with a stiffness much lower than 
predicted by the relevant Gibson-Ashby relation. 

• There is a pronounced stiffening in the initial stages 
of compression. 


It is also remarkable that both, the effective macro¬ 
scopic stiffness and the flow stress, agree closely between 
experiment and simulation, even though the ligament size 
of the experimental sample is tenfold larger than that of 
the simulation. This difference highlights that, in spite 
of several remarkable similarities between the deforma¬ 
tion phenomena in simulation and experiment, it would 
be inappropriate to over-emphasize the match between the 
results from the two approaches. Interatomic potentials 
in MD simulation do not simultaneously provide precise 
agreement with the surface tensions, surface stresses, and 
stacking fault energies of the experiment, nor do the ac¬ 
cessible timescales of MD allow for a reproduction of the 
orders of magnitude slower processes in the experiment. 
Nonetheless, discussing the simulation results on their own 
in a self-consistent way affords insights into the micro¬ 
scopic processes and their interactions by which nanoporous 
solids with microstructure similar - though not identical 
- to the real world deform. Based on a critical discus¬ 
sion of these processes in relation to the experiment, our 
study arrives at the following suggestions with respect to 
the acting deformation mechanisms: 

• Even in the absence of external load there is local 
yielding driven by capillary forces. As-prepared sam¬ 
ples therefore already exhibit lattice defects. 

• We propose that the surface-induced stresses create 
a distribution of shear stress that, in its tail, reaches 
the theoretical shear strength. The regions of some¬ 
what lesser shear stress are close to an elastic insta¬ 
bility, and are distinguished by vanishing shear stiff¬ 
ness. This is a key contribution to the anomalous 
compliance of NPG. 

• Dislocation activity occurs immediately upon load¬ 
ing, and compressive plastic deformation leads to the 
continuous accumulation of lattice dislocations. 

• The storage of lattice defects and our observations of 
locks imply that, in addition to densification, stan¬ 
dard dislocation-based mechanisms must contribute 
to the strengthening of NPG under compression. The 
dislocation-starvation scenario, which is often con¬ 
sidered in models for nano-object deformation, does 
not apply to NPG. 

• Similarly, the evolution of effective Young’s modulus 
during compression is not a simple density-related ef¬ 
fect. Other factors, including the relaxation of highly 
stressed (and, therefore, nearly shear-unstable) re¬ 
gions by faulting, are important. 
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Appendix A. Supplementary data 

Supplementary data associated with this article can 
be found in the online version. For the sake of clarity, 
the animation Movie-3 only shows reconstructed surface, 
planar defects (red), and dislocations (green). 
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